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Abstract 

We consider the impact of varying choices (and scales) on each side of the so-called "match- 
ing scale" in MLM-matched matrix-element + parton-shower predictions of collider observables. We 
explain how inconsistent prescriptions can lead to counter-intuitive results and present a few explicit 
examples, focusing mostly on W/Z+ jets processes. We give a specific prescription for how to improve 
the consistency of the matching and also address how to perform consistent tune variations (e.g., of the 
renormalization scale) around a central choice. Comparisons to several collider processes are included 
to illustrate the properties of the resulting improved matching, relying on AlpGen + Pythia 6, with the 
latter using the so-called Perugia 201 1 times, developed as part of this effort. 



1 Introduction 



The theoretical description of muhijet production in hadronic collisions is one of the key ingredients for 
the interpretation of the data from high-energy hadron colliders, the 1.96 TeV proton- antiproton Tevatron 
collider at Fermilab, and the proton-proton Large Hadron Collider (LHC) at CERN. Final states with multi- 
jets, possibly associated with electroweak gauge bosons, are in fact the dominant signature of the decay of 
heavy particles produced at high energy, whether in the Standard Model (top quarks and Higgs bosons), or 
in theories beyond the Standard Model (BSM), such as supersymmetry. The identification of these particles, 
and the study of their properties, requires an accurate modelling of the Standard Model (SM) sources of 
multijets. Great progress was achieved towards this goal in the past decade. On one side, the calculation 
of inclusive, parton-level, cross-sections to next-to-leading-order (NLO) in QCD has produced results for 
processes as complex as W-i-4 jets |[T|. On the other, algorithms have been developed and implemented in 
numerical codes to provide a complete description of the hadronic final states emerging from processes with 
up to 6 jets, merging the exact leading-order (LO) calculation of the partonic matrix elements (ME) with 
the evolution, provided by so-called shower Monte Carlo (MC) codes, of the partonic shower (PS) and the 
subsequent hadronization of the partons in to physical hadrons. 

The development of theoretical tools has been accompanied by experimental measurements, which pro- 
vide the necessary vahdation test-bed for these calculations. Parton-level NLO calculations provide a first- 
principle description of inclusive final states: they have an intrinsic high degree of precision, due to the 
reduced dependence on the unphysical choice of a renormalization and factorization scale and, furthermore, 
are not subject to modelling uncertainties related to the details of the non-perturbative phase of the final 
state evolution. Calculations based on the merging of LO matrix elements, shower evolution and hadroniza- 
tion, on the other hand, while affected by the larger scale-setting uncertainty due to the LO approximation, 
provide a fully exclusive description of the final states, and are therefore more suitable for the experimental 
analyses. Their ultimate goal is not only to give reliable estimates of the inclusive jet rates and energy distri- 
butions, but also to reproduce properties of the final states such as the jet inner structure and the distribution 
of softer particles produced outside of the jets, including those resulting from the evolution of the fragments 
of the original colliding hadron^ These properties, which depend on the details of the non-perturbative 
dynamics, can only be described through the phenomenological models embedded in the shower MC codes. 
The parameters of these phenomenological models need to be tuned using experimental data of some suit- 
able observables. The factorization assumption built into any description of large-Q processes justifies the 
use of these same parameters in the prediction of different observables, and provides the basis for the pre- 
dictive power of such tools. This assumption however must be validated with a direct comparison with data. 
Elements that need to be probed include the scaling with beam energy of the UE parameters, the universality 
of the parameters controlling the shower evolution and hadronization, and the overall independence of all 
parameters on the type of hard process. Deviations from the expected universality would highlight faults in 
the underlying modelling of effects beyond perturbative physics, or could be due to the insufficient precision 
of the perturbative description, in case NLO effects were to modify significantly the LO predictions. Differ- 
ences compatible with the theoretical systematics of the LO approximation could however be reabsorbed by 
modifying the perturbative parameters that govern the LO systematics, for example the renormalization and 
factorization scales, or the matching variables used in the matrix-element/shower merging algorithm. 

It is therefore important to understand the correlations between the effects of changing the soft and UE 
parameters on one side, and the perturbative parameters on the other. In this paper we present studies which 

'We refer to the ensemble of these particles as the underlying event, or UE 
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demonstrate that, in the tuning of ME-PS matched predictions, it is vital that there is consistency in the 
treatment of as in both the ME and PS components. While this is a general issue for all shower MCs, we 
consider as an expUcit example the merging of LO matrix elements with the Pythia 6 shower MC 0, as 
implemented in the framework of the AlpGen code |[3|, one of the reference tools for experimental multijet 
studies at the Tevatron and at the LHC. The most recent versions of Pythia (6.425) and AlpGen (2.14) codes 
were used for producing the results. 

On the Pythia 6 side, we consider several different tune variations of the interleaved ;?T-ordered parton- 
shower model ^4j|, focusing on the so-called "Perugia" set of tunes of |[5j|6|, ranging from the Perugia tune 
(from 2009) to the Perugia 20 11 updates that have been developed as part of this work, including systematic 
up/down variations of the shower activity (see the Appendix and |[5][6| for details). We also compare to the 
"DW" tune |7| of the virtuality-ordered shower model 18", '91. For Herwig fTOl, we include the "Jimmy" 
underlying-event model [11], with default parameters. We emphasize that the qualitative conclusions pre- 
sented in this paper carry over to other shower models, including the ones implemented in Pythia 8 |T2] 13 1 



and Herwig |14|, but the quantitative aspects should still be considered limited to the particular tunes 
and shower models studied here. We rely on Fastjet 1 15 ] for jet clustering and have further used the Rivet- 
based [ 16] mcplots web site [ 17] for some of our comparisons. 

The paper is organised as follows: in Section |2] we describe in detail the theoretical nature of the as consis- 
tency problem, and give a practical example of how it can be manifest in the prediction of high pj observ- 
ables. In Section|3]we show how a simple prescription can be applied to stabilise ME-PS tunings against this 
problem, propose a new tune for AlpGen -i- Pythia 6 matched predictions, and demonstrate the behaviour 
of this tune under tuning variations. In Section |4] we show that this new AlpGen -i- Pythia 6 tune is able 
to reproduce (within statistical errors) the Tevatron and LHC vector boson plus jets data. In addition we 
also the tune predictions to the jet shape measurements at the Tevatron and LHC. Finally we conclude in 
Section|5l 



2 The Importance of Consistent 0:5 Treatment in ME-PS Matched Predic- 
tions 

In this section we demonstrate that consistent treatment of as in ME-PS matched predictions is important 
in order to achieve the desired accuracy in the prediction of high pj observables. We first present the 
theoretical arguments behind this, and then go on to show and explain that without adopting this approach 
one can observe undesirable and counter-intuitive effects on experimental observables. 



2.1 Theoretical Background 

The philosophy behind matching prescriptions such as the MLM one p"8p9} employed by AlpGen is to sep- 
arate phase space cleanly into two distinct regions; a short-distance one, which is supposed to be described 
by matrix elements, and a long-distance one described by parton showers. In the long-distance region, real 
and virtual corrections, with the latter represented by Sudakov factors, are both generated by the shower and 



are intimately related by unitarity (for pedagogical reviews, see, e.g., |20 21 1). On the short-distance side, 
the real corrections aie generated by the matrix elements while the virtual ones are still generated by the 
shower. 
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Much effort has gone into ensuring that the behaviour across the boundary between the two regions be as 
smooth as possible. CKKW showed | |22J that it is possible to remove any dependence on this "matching 
scale" at NLL precision by careful choices of all ingredients in the matching; technical details of the im- 
plementation are important, and the dependence on the unphysical matching scale may be larger than NLL 
unless the implementation matches the theoretical algorithm precisely 



Especially when two different computer codes are used for matrix elements and showering, respectively (as 
when AlpGen or MadGraph |26 | is combined with Pythia 6 or Herwig), inconsistent parameter sets between 
the two codes can jeopardise the consistency of the calculation and lead to unexpected results, as will be 
illustrated in the following sections. 

To give a very simple theoretical example, suppose a matched matrix-element generator (MG) uses a differ- 
ent definition of oc^ than the parton-shower generator (SG). Suppressing parton luminosity factors to avoid 
clutter, the real corrections, integrated over the hard part of phase space, for some arbitrary final state F, will 
then have the form 

4ti= TdOf+iaf^^lMf+ip, (1) 

where we have factored out the coupling corresponding to the "-i-l" parton and suppressed the dependence 
on any other couplings that may be present in \Mp^i p. The virtual corrections at the same order, generated 
by the shower off F, will have the form 

^exci ^ ^mci _ J ^ dz £ ^/^.(z) \Mp\' + ff{a^) , (2) 

with Pi{z) the DGLAP splitting kernels (or equivalent radiation functions in dipole or antenna shower ap- 
proaches). If the two codes use the same definitions for the strong coupling, af^ = a]^*-", then the fact 
that P{z)/Q^ captures the leading singularities of |M/r+ip guarantees that the difference between the two 
expressions can at most be a non-singular term. Integrated over phase space, such a term merely leads to 
a finite ^(ots) change to the total cross section, which is within the expected precision. Indeed, it is a cen- 
tral ingredient in both the MLM and (L)-CKKW matching prescriptions that a reweighting of the matched 
matrix elements be performed in order to ensure that the scales appearing in match smoothly between 
the hard and soft regions. Thus, we may assume that the choice of renormalization scale after matching is 
}i ^ Pt on both sides of the matching scale, where is a scale characterising the momentum transfer at 
each emission vertex, as established by |27]28j and encoded in the CKKW formalism |22 1. 



In the case of the CKKW approach as implemented in the Sherpa MC framework p9| , this prescription can 
be controlled exactly, since the matrix element and the shower evolution are part of the same computer code 
and hence naturally use the same definition. This is also true in Lonnblad's variant [ )23J of the algorithm, 
used in Ariadne [[301. In the case of codes like AlpGen or Madgraph, on the other hand, an issue emerges. 
These codes are designed to generate parton-level event samples to be used with an arbitrary shower MC. 
Different shower MCs however use slightly different scales for the paiton branchings, as a result of different 
approaches to the shower evolution, and may use different values of Aqcd> as a result of the tuning of the 
showers and/or underlying events. A possible mismatch therefore arises in the values of Us used by the 
matrix-element calculation and those used by the shower. 

If there is a mismatch in Aqcd or as{Mz), then this will effectively generate a real-virtual difference whose 
leading singularities are proportional to 

a^b,\n{^^y-^l^m\Mp\\ (3) 
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Table 1: The three cases, A, B, and C discussed in the text, for an arbitrary reference A value. 

which is of next-to-leading logarithmic order (unless Amg ~ ^sg, in which case it vanishes). Similarly, even 
if both matrix-element and shower codes are using the same Aqcd, but they use different running orders, 
then there will be an iff {a^ \n{p^ / A^)) mismatch, which may also become large if p-j- » A. 

To be more concrete, let us consider a specific example. Compare A) a matched MG-i-SG calculation which 
uses the same Aqcd value on both sides of the matching to B) a calculation in which the value used on the 
MG side is reduced to half its previous value but the SG one remains the same, as summarised by the two 
first columns of tab.[T] Going from case A to B, the following changes result: 

1 . The number of (F + 1 ) states added by the MG decreases, due to the lowering of the Aqcd value on 
the MG side, while the number of surviving F states remains constant, since the shower Sudakov is 
not modified. The total estimated cross section therefore decreases. 

2. At the differential level, the smaller number of (F + 1) states combined with the unchanging number 
of F states implies smaller absolute jet cross sections and smaller fractions ajet/cJtot. 

Similarly we may consider what happens if C) we reduce the Aqcd value on the SG side instead, as sum- 
marised in the last column of tab.[T] Going from case A to C, the following changes result: 

1 . The number of (F + 1 ) states added by the MG remains constant, while the number of surviving 
F states increases, since the SG is generating fewer branchings. The total estimated cross section 
therefore increases. 

2. Since the number of (F + 1) states is constant, while the shower is made less active, the final jets will 
actually be narrower, which increases the rate of reconstructed jets at any given fixed pj value. 

3. Since both the total cross section increases and the number of reconstructed additional jets also in- 
creases, jet fractions can either increase or decrease. 

In particular, note the somewhat counter-intuitive effect that decreasing the shower oc^ value actually in- 
creases the jet rates in a matched calculation, while it normally decreases them in a standalone shower 
calculation. 

Since, as was discussed above, inconsistencies among the choices on the two sides can lead to differences at 
the NLL level, it is obviously important to ensure that they are consistent within a reasonable margin. This 
is particularly true in the context of event-generator tuning, in which specifically the NLL components of 
the shower description are sought to be optimized with respect to measured data, and hence changes at this 
level could effectively destroy the tuning. 

Finally, we remind the reader that a change in Aqcd can be interpreted as a change in the opposite direction 
of the renormalization scale argument (for constant Aqcd)> modulo small flavour threshold effects that we 



4 



5 


1.8 


Q 




■»-' 


1 .6 


□_ 
+ 


-H A 

1 .4 


Alp. 


1 .2 


o 


1 


o 




To 


0.8 


DC 






0.6 




0.4 



p p;\''s=1 .96 TeV; (W Gv)4^jets 



tn^|<1 .1 ; p^>20 GeV; ^^>30 GeV \ 



4= 



Pythiastat. err. 
Alp.+Her. /Alp.+Pyt. DW 
Alp -i-Her.Jim. .''Alp.4Pyt. DW 
Alp.+Pyt. PO/Alp.-i-Pyl. DW 



30 40 50 



10^ 2x10^ 

jet [GeV] 



Figure 1 : Ratio of predictions for the leading-jet Ej spectrum in W+jets final states at the Tevatron, 
obtained with AlpGen plus various MC codes and tunes. The leading jet observable is defined at 
the particle-level as in the CDF W-i-jets analysis |[32| . 



shall ignore here. This is easy to realise from the definition of the coupling, 



2< 1— loop 



1 



1 



(4) 



Thus, we may write renormalization scale variations (e.g., by a factor of 2 in each direction) either by 
applying a prefactor directly on the renormalization scale argument of or by applying the inverse of that 
factor to Aqcd while keeping the renormalization scale argument unchanged. Due to the technical structure 
of the codes, the former is more convenient in AlpGen (via the ktf ac setting) whilst the latter is more 
convenient for Pythia 6. 



2.2 Examples of the interplay between tunes and matching 

In this section we give several examples of how the issues in ME-PS matching described in Section |2] can 
affect high-jOx observables using AlpGen interfaced to Pythia 6 with DW Perugia (PO) ||5j|3Tj and 
Perugia 2010 (P2010) tunes ||6). 

In Fig. [T] we show the ratio of predictions for the transverse energy (Et) spectrum of the leading jet (that 
jet with the highest Ej per event) in W-njet final states at the Tevatron, obtained by the merging of AlpGen 
with different shower codes; Herwig, Pythia 6 virtuality-ordered shower (DW), and /^T-ordered shower 
(Perugia 0). The differences between Herwig and Herwig plus Jimmy at small Ej can be explained by 
the different amounts of energy that, in the various cases, are deposited by the UE in the jet cones. In 
particular, as shown in Fig.[2j these differences can accommodate the slight shape discrepancy between data 
and AlpGen -i-Herwig that was noted, at small Ej, in the CDF study [32]. It is difficult, however, to attribute 
to the UE energy the significant differences seen in Fig. [T] at large Ej. 
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Figure 2: Comparison of CDF data |32| with the leading-jet Ex spectrum predicted by AlpGen 
plus various MC codes and tunes. 



In order to investigate the source of the differences in the predictions, systematic parameter variations of the 
perturbative and non-perturbative model components of Pythia 6 have been studied using the Perugia family 
of Pythia 6 tunes with Perugia as the central tune and Perugia Hard and Perugia Soft as the systematic 
variation tunes. The Perugia Soft and Perugia Hard tunes both use the same Parton Density Function (PDF) 
as Perugia 0, CTEQ5L j33| , but differ in the values of Pythia 6 parameters controlling both perturbative and 
non-perturbative activity levels. In comparison to Perugia Soft, the Perugia Hard tune has more perturbative 
(initial and final state radiation) activity but less non-perturbative (multiple interactions, beam remnant and 
hadronization) activity. Perugia Soft on the other hand has less perturbative but more non-perturbative 
activity than the Perugia tune. In order to investigate the interplay of the tuning variations with the MLM 
matching, the effect of the change of the tune on the physics observables in both the Pythia 6 standalone 
case and AlpGen -i- Pythia 6 case is presented. 

In Fig.|3]the distribution of jet multiplicity (Njet) in W-f-jets events is compared for events generated with the 
Perugia 0, Perugia Hard and Perugia Soft tunes. The Njet observable is defined at the particle-level according 
to the definition used in the ATLAS measurement of the W-f-jets cross-section at ^/s=l TeV [34] . Jets are 
clustered from stable particles using the anti-Kt jet algorithm [35] with the radius parameter R = 0.4, and 
considered in case they satisfy the following kinematic cuts: pj > 20 GeV and |t]| < 2.8. Comparisons 
are performed for both the Pythia 6 standalone (left) and AlpGen -i- Pythia 6 (right) cases. For the Pythia 6 
standalone case we observe that Perugia Hard tune yields more high-jOx jets than the Perugia tune while 
Perugia Soft yields less final state jets correspondingly. For the AlpGen -i- Pythia 6 case an opposite trend is 
observed: Perugia Hard tune yields less high-7?T jets and Perugia Soft tune yields more high-px jets. In order 
to determine which modelling components of the Perugia Soft and Perugia Hard tunes cause this behaviour, 
we considered the effect of varying individual sets of parameters of the Perugia Soft and Perugia Hard tunes 
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in AlpGen + Pythia 6 predictions. Parameters were grouped according to the modelling aspect they control 
into Initial State Radiation (ISR), Final State Radiation including the FSR from the ISR partons (FISR), the 
Underlying Event (UE) and Colour Reconnections (CR) block^ 

Dedicated samples where only parameters of an individual block were varied in the ranges used in Peru- 
gia Soft and Perugia Hard tunes on top of the Perugia tune were produced. The results of the study, in 
terms of the cross-section contribution of each AlpGen sub-sample (after MLM matching), are given in 
Table |2] As was already noted in Fig. [3j the cross-section for multijet production in the Perugia Hard case 
decreases with respect to Perugia 0, and vice versa for Perugia Soft. From Table |2] we see that the param- 
eter blocks that produce this affect are the ISR and FISR blocks, while the impact of the CR and UE block 
variations on the cross-sections is negligible. In addition to the simultaneous variations of parameters in 
the blocks, we have also performed individual parameter variations for each of the parameters in order to 
check that potential correlations between the parameters do not affect the conclusions. Studies have also 
been performed for the Hadronization and Beam Remnant blocks of 16 ]. The variations of these parameters 
also had a negligible effect on the kinematic distributions and cross-section values. 

In Fig. |4] we demonstrate that the increased parton shower activity can indeed lead to the reduced cross- 
section (and softer jet spectra) due to the increased rates at which the AlpGen -i- Pythia 6 events are vetoed 
during the MLM matching. In the figure the distributions of the events that pass or fail (ISVETO=0 or 
ISVETOt^O) the MLM matching criterion are shown for the exclusive sub-sample of AlpGen -i- Pythia 6 
Perugia 2010 W-i-jets events with exactly three additional partons from the matrix element in the final state|^ 
Each of the distributions is normalised to unit area. The distributions are shown as a function of the largest 
Pi shower emission from the initial state radiation (left) and as a function of the largest multiple proton- 
proton interaction. In the left hand side figure we see that the events are rejected with higher probability, the 
larger the p^ of the hardest ISR branching in the event. Therefore, a Pythia 6 standalone tune which increases 
the ISR activity can, somewhat counter-intuitively, reduce the rate for multijet and hard emissions. In the 
right hand side we demonstrate that the events are accepted and rejected independently of the transverse 
momentum of the hardest multiple interaction in the event (which is the desired behaviour of the matching 
application used with the parton shower code). 

To conclude, the origin of the differences observed in the predictions of tunes with different ISR/FSR activity 
matched to AlpGen is rather due to the mismatch between the jet-emission probability predicted by the 
matrix elements and by the shower. This comes from the mismatch in the value of as discussed earlier, 
arising from different values of Aqcd or from the use of a different evolution variable in the shower. If the 
value of as in the shower increases, the emission rate of additional jets during the shower evolution will 
increase. Since the matching algorithm rejects events with extra jets generated by the shower, to replace 
them with events where the jet is accounted for by a higher-order matrix element calculation, a larger value 
of as in the shower leads to a higher rejection rate. Unless this change in as is accompanied by a similar- 
change in the matrix element calculation, the additional rejection is not compensated by the relative increase 
in rate for the higher-order parton-level contributions, leading to the effects reported in this section. 

This important interplay between MC parameters, which are typically tuned to "soft" observables such as 
UE or the small-/7T DY spectrum, and the performance of the matching algorithms for "hard" observables, 
calls for particular attention when adopting new UE tunes in the framework of multijet studies with matrix- 
element matching. Along the same lines, it should be kept in mind that, tuning a stand-alone shower MC to 



A.2 



The parameter blocks organisation is similar to the one introduced in |j6| and are listed in 
^The observations in the text are largely independent on the final state parton multiplicity 
"^These pj values are reported by Pythia 6 parameters VINT (357) (ISR) and VINT (359) (MPI) respectively. 
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Figure 3: Jet {p^ > 20 GeV and < 2.8) multiplicity distribution in W+jets electron channel 
events in pp collisions at 7 TeV. Distributions are shown for the samples generated with Pythia 6 
standalone (left) and with AlpGen + Pythia 6 (right). For each case the distributions obtained 
when using Perugia (PO), Perugia Hard (Phard) and Perugia Soft (Psoft) tunes are shown. All 
distributions are scaled so that the value of the first bin agrees with the ATLAS measurement p4) . 



tune 


NpO 


Npl 


Np2 


Np3 


Np4+ 


total [pb] 


Phard 


7287 ± 3.9 


728 ± 2.6 


141 ± 1.3 


27 ± 0.2 


6.6 ± 0.2 


8190 ± 8 


PO 


7556± 3.6 


814 ±2.7 


166 ± 1.3 


32 ± 0.3 


7.8 ±0.3 


8576 ± 8 


Psoft 


7804 ± 3.4 


944 ± 2.8 


207 ± 1.5 


42 ± 0.3 


10.1± 0.3 


9007 ± 8 


PO with Phard ISR 


7207 ± 6.9 


735 ± 2.6 


143 ± 1.3 


27 ± 0.2 


6.9 ± 0.2 


8119 ± 11 


PO with Psoft ISR 


7831 ±4.9 


881 ± 2.7 


186 ± 1.4 


36 ± 0.3 


8.8 ±0.3 


8943 ± 10 


PO with Phard FISR 


7548 ± 6.0 


814 ±2.7 


167 ± 1.3 


32 ± 0.3 


7.8 ±0.3 


8569 ± 10 


PO with Psoft FISR 


7505 ± 6.1 


878 ± 2.7 


188 ± 1.4 


37 ± 0.3 


9.4 ± 0.3 


8617 ± 10 


PO with Phard UE 


7513 ± 6.1 


826 ± 2.7 


171 ± 1.4 


33 ± 0.3 


7.8 ±0.3 


8551 ± 10 


PO with Psoft UE 


7576 ± 5.9 


817± 27 


166 ± 1.3 


32 ± 0.3 


8.1 ±0.3 


8599 ± 10 


PO with Phard CR 


7561 ± 5.9 


821 ± 2.7 


167 ± 1.3 


32 ± 0.3 


8.1 ±0.3 


8589 ± 10 


PO with Psoft CR 


7556 ± 5.9 


815 ± 2.7 


165 ± 1.3 


32 ± 0.3 


8.1 ±0.3 


8576 ± 10 



Table 2: Impact of different variations of the Pythia 6 Perugia 2010 tunes on the cross sections of Alp- 
Gen W+jets sub-samples with different matrix element parton multiplicities, and the total inclusive W cross 
section. For the studies sub-samples with up to four additional partons from the matrix element were gener- 
ated. The matching is performed inclusively for the highest parton multiplicity sub-sample and exclusively 
for other sub-samples. The tabulated cross-sections were extracted after the MLM matching and parton 
shower. The errors shown are statistical only. The parameter settings of the various setups are discussed in 
the text. 
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Figure 4: Distribution of the probabilities for the event acceptance (ISVETO=0) or rejection 
(IS VETO/0) during the MLM matching step, as a function of the largest p-^ shower emission from 
the initial state radiation (left) and the largest multiple proton-proton interaction in the event 
(right). The events were generated using exclusive sub-sample of AlpGen -i- Pythia 6 Perugia 2010 
VK-f-jets events with exactly three additional partons from the matrix element in the final state for 
pp collisions at 7 TeV and AlpGen -i- Pythia 6 Perugia 2010 tune. 



better model multijet final states, will force it to emulate effects present in the multiparton matrix elements. 
Using such tunes with matrix-element matching therefore requires ad-hoc modifications of the matching 
algorithm, or of its parameters. 



3 Stabilising ME-PS Matched Timings 

In this section we discuss how to overcome the problems discussed in the previous section with a simple 
prescription, and outline a tuning strategy that should allow to consistently optimize, in the context of the 
Pythia 6 shower MC, the description of both the UE and the high-Ej properties of final states. 



3.1 A New AlpGen + Pythia 6 as Consistent Tune 



As it was explained in Section 2.1 and practically demonstrated in Section [2^2} it is highly desirable to have 



a consistent treatment of as on either side of the ME and PS boundary. In Appendix A.l| the relevant settings 



for a new as consistent AlpGen -i- Pythia 6 tune are described in detail. In this tune, the as consistency is 
essentially ensured by setting the effective value of Aqcd to be the same throughout the Pythia 6 parton 
shower algorithms and in the AlpGen matrix elements. A consistent choice for Aqcd of 

Ag,^D~0.26, (5) 

is made, where the superscript indicates the number of flavours. This choice is informed by comprehensive 
Professor tunings | [36l[37| of the p^-ordered shower in Pythia 6 Q to event shapes and other LEP data. Note 
that the settings for Pythia 6 are those of the central Perugia 201 1 (P201 1) tune ||6|, which was inspired by 
these studies. We will refer to this new tune of AlpGen -i- Pythia 6 as the Perugia 201 1 "matched" tune. 
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Figure 5: Comparison of AlpGen + Pythia 6 (pj >20 GeV) jet multiplicity (left) and leading jet 
transverse momentum (right) distributions in W+jets electron channel events. The samples are 
generated using different AlpGen + Pythia 6 parameter setups described in the text. 



3.2 Tests of the Consistent as Approach: Behaviour Under Scale Variations 



In this section we study the behaviour of the new AlpGen + Pythia 6 Perugia 2011 "matched" tune under 
Aqcd variations to demonstrate that, with a consistent treatment of as, the expected behaviour of ME-PS 
matched predictions under variations of tuning parameters is restored. VK+jets events selected with the same 
criteria applied for Fig. [3] are used. Figure [5] shows the jet multiplicity (left) and leading jet transverse mo- 
mentum (right) distributions for the Perugia 201 1 "matched" tune and four variant tune samples generated 
with different Aqcd values. Two samples, labelled as "A Alp. t" and "A Alp. J,", have Aqcd respectively 
increased and decreased by a factor of 2 only in the ME calculation. This is achieved by setting respectively 
the AlpGen parameter ktfac to 1/2 and 2. The increase (decrease) of the Aqcd value in AlpGen results 
in more (less) jets and a harder (softer) leading jet spectrum as shown in Fig. [5] The two samples labelled as 
"A PS ti ^ Alp. t" and "A PS |, A Alp. |" correspond to a consistent variation of Aqcd both in the ME 
and PS, with Aqcd respectively increased and decreased by a factor of 2. The impact of these variations is 
qualitatively similar to the case where Aqcd is only varied in the ME, restoring the expected behaviour of 
ME-PS matched prediction under variation of Aqcd- However, the samples with Aqcd varied simultane- 
ously in the ME and in the PS exhibit a smaller deviation from the nominal sample. The mitigation of the 
impact of a Aqcd coherent change in a ME-PS matched sample compared to the same change only in the 
ME calculation is due to the interplay between the radiation produced by PS and the matching algorithm, as 
detailed in Section 2. 1 While the choice of the xlclu parameter allows to directly adapt AlpGen to possible 
future changes in the choice of Aqcd in Pythia, the variation of the ktfac parameter in the standard range 
0.5 < ktfac < 2 can be used to estabhsh the range of the systematical uncertainty, or to tune the description 
of specific observables. 
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4 Comparisons with Data 



In this section we demonstrate that the new AQco-consistent Perugia 201 1 "matched" tuning of AlpGen + Pythia 6 
introduced in Section 3.1 compares well with recent Tevatron and LHC measurements, and that, with the 
arrival of improved precision measurements, there should be room for further tuning of these predictions. 



4.1 Z/W+ jets production 

The figures that follow show comparisons of AlpGen + Pythia 6 Monte Carlo predictions to measurements of 
pubhshed Z+jets and W+jets processes from CDF | [32l[38|[39) and W+jets from ATLAS [40f| These cross- 
section measurements are corrected for all known detector effects to particle level and compared to Monte 
Carlo predictions. The Monte Carlo predictions of V+jet production cross-sections are formed by clustering 
the stable final state particles (t > 10 ps) following parton shower and hadronization of the unweighted 
events. This clustering is done using the same jet algorithm as the measurement, as implemented in the 



Fastjet 1 15 1 package. All stable final state particles are used, with the exception of the leptons that result 
from the decay of the signal W or Z boson. The decay leptons are corrected for the final state QED radiation 
such that their 4-momentum is equivalent to that before radiation. After the events have been clustered, the 
restrictions on the allowed phase space of the jets and of the W/Z boson decay products are applied to be 
consistent with the measurement to which we are comparing. The prediction of the final AlpGen -i- Pythia 6 
cross-sections contains contributions from V-i-O, 1, 2, 3 parton samples (showered with exclusive MLM 
matching), and V+4 parton samples (showered with inclusive MLM matching). 

Two different AlpGen -i- Pythia 6 generations are compared to the data; the new AlpGen -i- Pythia 6 Peru- 
gia 201 1 "matched" tune introduced in Section[3?T](labelled "Alp.-i-Pyt. P201 1"), and an AlpGen -i- Pythia 6 
prediction using the default settings of AlpGen and the Pythia 6 DW tune (labelled "Alp.-i-Pyt. DW"). The 
ratio of the matched predictions to the data are shown and compared. Additionally, the results of variations 
of ktf ac by factors of 0.5 and 2.0 in the AlpGen -i- Pythia 6 Perugia 2011 "matched" prediction are shown 
as solid lines. The hatched regions show the total error (statistical plus systematic) propagated to the the- 
ory/data ratio from the data measurements. The error bars on the points show the statistical error on the 
theoretical prediction. 

In Figure|6]we show the ratio of the predicted theory cross-sections to the data for the CDF Z-f-jets measure- 
ment | [39} . In this measurement jets are defined by the CDF midpoint algorithm (46], with Rcone = 0.7 and 
are required to have p^^' > 30 GeV and \y^^'\ < 2.1. In Figure]? we show the ratio of the predicted theory 
cross-sections to the data for the CDF W-i-jets measurement |32 . In this measurement jets are defined by 
the CDF JetClu algorithm |47|, with Rame = 0.4 and are required to have Pt > 20 GeV and |t]| < 2.5. In 



Figure we show the ratio of the predicted theory cross-sections to the data for the ATLAS W-f-jets mea- 



surements [40]. In this measurement jets are defined by the anti-Kt algorithm |35 1, with a radius parameter 
R = 0.4 and are required to have pj > 20 GeV and | T] | < 2.8. 

The AlpGen -i- Pythia 6 Perugia 201 1 "matched" prediction compares well with the measured cross sections 
both as a function of the inclusive jet multiplicity and jet pj. In particular, the prediction correctly describes 
the low Pt- region of the differential cross section and the jet sub-structure without presenting any significant 
disagreement with data at high pj. This shows that it is possible to tune separately the long- and short- 



^ Measurements of these processes at the Tevatron have also been performed by DO |4llj43|. Preliminary, higher statistics, 
studies of IV /Z+jets production at the LHC have been reported by ATLAS (34|[44) , and CMS |45| . 
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Figure 6: (Top) The ratio of predicted theory and CDF measured data cross-sections for the 
production of a Z — )• boson in association with at least Njet jets |[39j. In the left hand figure the 
theory predictions are not normahsed to the data. In the right hand figure the theory predictions 
are normalised such that they equal the data measurement in the > 1 jet bin. (Bottom) The ratio of 
predicted theory and CDF measured data cross-sections for the production of a Z — boson in 
association with at least 1 jets (left hand side) and at least 2 jets (right hand side) as a function of 
jet Ej. In the left hand plot, the theory prediction is normalised such that the predicted rate for > 1 
jet production is equal to that measured in the data. In the right hand plot, the theory prediction is 
normalised such that the predicted rate for > 2 jet production is equal to that measured in the data. 



distance contributions of the prediction to obtain a satisfactory description of the observables in the whole 
experimental accessible phase space. Remarkably, the prediction describes the data both at ^fs = 1 .96 TeV 
and y/s = 1 TeV. This illustrates that the scaling properties of the long-distance contribution with the centre 
of mass energy of the collision does not produce unexpected effects in the high region of the cross section. 
A coherent rescaling of as with ktf ac =0.5, 2 has a little effect on the shapes of the differential cross- 
sections, while for the inclusive Njet cross-sections it produces variations that bracket the default prediction. 
The ktf ac parameter can therefore be used to explore the sensitivity of the prediction to a variation of the 
renormalization and factorization scale, other than allow tuning on data. With the statistics of the currently 
available measurements there is no much room for optimising the parameters of the AlpGen -i- Pythia 6 
Perugia 2011 "matched" prediction in a tune that better describes the measurements. However, with the 
imminent arrival of more precise higher-statistics LHC measurements, this should be possible in the near 
future. 

Even though not explicitly shown here, in the relevant publications one can find similarly good agreement 



between the available measurements ||32}[34||44J and predictions based on AlpGen -i- Herwig. 
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Figure 7: (Top) The ratio of predicted theory and CDF measured data cross-sections for the pro- 
duction of a W — ev boson in association with at least Njet jets 1 32 1 . In the left hand figure the 
theory predictions are not normalised to the data. In the right hand figure the theory predictions 
are normalised such that they equal the data measurement in the > 1 jet bin. (Bottom) The ratio 
of predicted theory and CDF measured data cross-sections for the production of events containing 
a W — ev boson in association with at least 1 jets (left hand side) and at least 2 jets (right hand 
side), as a function of the leading jet Ej (left hand side), and the sub-leading jet Et (right hand 
side). In the left hand plot, the theory prediction is normalised such that the predicted rate for > 1 
jet production is equal to that measured in the data. In the right hand plot, the theory prediction is 
normalised such that the predicted rate for > 2 jet production is equal to that measured in the data. 
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Figure 8: (Top) The ratio of predicted theory and ATLAS pubUshed cross-sections |40| for the 
production of a W — boson in association with at least Njet jets. In the left hand figure the 
theory predictions are not normalised to the data. In the right hand figure the theory predictions 
are normalised such that they equal the data measurement in the > jet bin. (Bottom) The ratio of 
predicted theory and ATLAS data cross-sections for the production of events containing a — )• eV 
boson in association with at least 1 jets (left hand side) and at least 2 jets (right hand side), as a 
function of the leading jet p-^ (left hand side), and the sub-leading jet p-^ (right hand side). In the 
left hand plot, the theory prediction is normalised such that the predicted rate for > 1 jet production 
is equal to that measured in the data. In the right hand plot, the theory prediction is normalised 
such that the predicted rate for > 2 jet production is equal to that measured in the data. 
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4.2 Jets shapes 



Finally we test the ability of the new AlpGen + Pythia 6 Perugia 2011 "matched" tune and the systematics 
variations to describe the jet shapes at the LHC and the Tevatron. 



For the LHC, the jet shapes measured in inclusive jet production by the ATLAS collaboration | [48 | are taken 
as reference. For this measurement the jets are reconstructed using the anti-kt algorithm with the distance 
parameter R=0.6, the transverse momentum range 30 GeV < p-j < 600 GeV and rapidity in the region \y\ < 
2.8. The jet shapes are expected to be sensitive to both perturbative (parton shower) and non-perturbative 
(fragmentation and underlying event) modelling aspects. We perform the data comparisons for both the 
AlpGen + Pythia 6 and Pythia 6 standalone cases. The samples are generated using different Pythia 6 
standalone and AlpGen + Pythia 6 parameter settings as follows: for Pythia 6 standalone the Perugia 2011 
and the associated systematics tunes Perugia 2011 radHi and Perugia 2011 radLo are compared. The same 
tunes are also used for the generating the AlpGen + Pythia 6 distributions whereby the Aqcd values are 
always set to the same values in AlpGen (using ktfac) and Pythia 6 (i.e. the Perugia 2011 "matched" 
central settings and the systematic variations around the central settings). The setups are compared to the 
integral jet shape distributions as measured in the data. The integral jet shape is defined as the average 
fraction of the jet px that lies inside a cone of radius r concentric with the jet cone pSj: 



xi>^r) = —y^^^,0<r<R. (6) 
The sum is performed over all the Njet jets in the kinematic region of interest. 

In Figure [9] the integral jet shape distributions are compared to the ATLAS data for the jets in the trans- 
verse momentum ranges of 40-60 GeV (top) and 260-310 GeV (bottom) in the whole measured rapidity 
range {\y\ < 2.8). We observe that both Pythia 6 standalone (left) and AlpGen -i- Pythia 6 (right) with Pe- 
rugia 201 1 provide reasonably good description of the jet shapes. Due to MLM matching the jets in the 
AlpGen -i- Pythia 6 case tend to be more narrow than in the Pythia 6 standalone case. 



For the Tevatron, the shapes of jets produced in association with a Z boson as measured by CDF [ |49| are 
used. In this measurement jets are defined by the CDF midpoint algorithm [46 ], with Rcone = 0.7 and are re- 



quired to have pj. > 30 GeV and \y^'^' \ < 2. 1 . Figure 10 shows good agreement between AlpGen -i- Pythia 6 



and the measurement for both the DW and Perugia 201 1 tunes. 



The comparisons in Figures |9] and 10 (as well as comparisons to the jet shapes in other kinematic regions 



and comparisons to further LHC measurements) reveal no major short-comings of the AlpGen -i- Pythia 6 
Perugia 2011 "matched" tune. The Perugia 2011 tune has been developed by tuning Pythia 6 standalone, 
whereby the effective value of Aqcd was set to be the same throughout the Pythia 6 parton shower in the 
anticipation of using it with the AlpGen matrix elements using the same effective Aqcd value. The agree- 
ment with the measured jet shapes data could therefore potentially be improved by performing a dedicated 
tuning of AlpGen -i- Pythia 6. 



5 Conclusions 



We have shown that, in the context of tuning ME-PS matched predictions, it is vital that the tuning adopted 
ensures a consistent treatment of as on either side of the "matching boundary". In the case of Alp- 
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Figure 9: Comparison of the integral jet shapes as measured by ATLAS [48] with the predic- 
tions of the Pythia 6 standalone (left) and AlpGen + Pythia 6 (right) using Perugia 2011, Peru- 
gia 2011 radHi and Perugia 2011 radLo tunes. The comparisons are performed for the jets with 
\y\ < 2.8 and pr ranges of 40-60 GeV (top) and 260-310 GeV (bottom). 
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Figure 10: A comparison of the predicted and CDF measured jet shape in Z+jets events, showing 
the fraction of the total energy of a jet of radius R contained within a radius r, ^{r), as a function 
of r/R. "Alp.+Pyt. DW" uses the default AlpGen parameters and the DW Pythia tune, "Alp.+Pyt. 
P20H" is described in Section 3.1 and it includes a consistent choice of Aqcd in AlpGen and 
Pythia. 



Gen + Pythia 6 matched predictions, we have outlined a simple prescription to ensure this. This can be 
easily generalised, and applied to the case of matching AlpGen with other shower MCs, such as Herwig. 
We have then given an example of such a tune that compares well to Tevatron and LHC measurements of 
vector boson plus multijet final states. In addition, we have shown how consistent variations around a central 
ME-PS matched tune can be performed, so as to define a systematic uncertainty on that prediction. This 
knowledge should prove valuable in defining a new set of consistent ME-PS tunes for the precise future 
study of LHC multijet final states. 
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A Appendix 



A.l An as Consistent AlpGen + Pythia 6 Tune 

We here describe the parameters in the AlpGen and Pythia 6 codes that are important for ensuring the 
consistency of matching, and give the settings of these parameters which describe a new as consistent 
AlpGen + Pythia 6 tune. Note that the settings for Pythia are those of the central Perugia-2011 tune [6|, 
which was inspired by these studies. 

A.1.1 AlpGen Parameters 

For AlpGen, the relevant parameters controlling the magnitude, renormalization scale, and running order 
for as are set to these values Q 

• The prefactor rescaling the value of the renormalization scale in the CKKW scale-setting procedure: 
ktfac = 1.0 

• Aqcd (5-flavour): xlclu =0.26 

• Running order: Ipclu = 1 

A.1.2 Pythia Parameters 

To the extent that the explanations given here are necessarily somewhat brief, we recommend the interested 
reader to follow up on the definitions of the parameters below in the program's comprehensive manual |j2|. 

The Strong Coupling in Pythia 6: there are several different ways of specifying the parameters control- 
ling as- In order for the user to have explicit control of all of them, we use the option MSTP(3)=1, which 
allows to specify formally independent Aqcd values for each of the different algorithmic components in 
Pythia. As a special case, we may then choose to set all those values equal, or at least set those equal that 
correspond to initial- and final-state radiation. These are: 

• PAR? (61) for Aqcd for ISR. 

• PARJ (81) for Aqcd for FSR inside resonance decays. 

• PARP(72) for Aqcd for FSR outside resonance decays (e.g., FSR off hard jets from the matrix element 
and/or from ISR). 

The number of flavours with which to interpret Aqcd must also be specified. Since a value of «j = 5 is 
hard-coded in at least one sub-algorithm of Pythia, we advise to always translate Aqcd values to 5-flavour 
ones and, correspondingly, set the Uf parameter MSTU(112)=5. 

^The parameters xlclu and Ipclu were introduced in connection with this work, and are implemented in AlpGen starting from 
v2.14. The default behaviour, if these parameters are not set, is to assign the values inherited from the PDF (as was the case for 
AlpGen versions before 2.14). 
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Further, in order to avoid that any of these values are modified by the code, we set MSTP(64)=2, and 
PARP (64) =1 . 0. The former forces the code to keep Aqcd unmodified for ISR. In particular, the translation 
from "MSbar" to "CMW", which is applied for MSTP(64)=3, is not performed. This is equivalent to inter- 
preting the effective Aqcd value as already being in a scheme similar to CMW. The latter, PARP (64) =1 . 0, 
sets the prefactor for the renormalization scale used for ISR equal to unity, i.e., the renormalization scale 
will just be Pj. Any re-interpretation of Aqcd> for instance to translate between different effective scheme 
definitions or to introduce multiplicative factors on the effective renormalization scale, for scale variation 
purposes, should then be imposed directly on the three Aqcd values above. This is the prescription followed 
in the so-called Perugia 2011 tunes which were developed as part of this effort, with parameters as listed 
int6J. 

Finally, one needs to settle on an effective value for Aqcd- According to comprehensive Professor tunings 
f36„37J of the px-ordered shower in Pythia Q to event shapes and other LEP data, one needs values of order 

AgD~0-26, (7) 

where the superscript indicates the number of flavours. We interpret this as an effective value, derived 
directly from data using a "Pythia scheme" that is defined numerically by Pythia's shower algorithm. It is 
not necessarily directly comparable to MS determination^ 



The Strong Coupling in Pythia 8: Although we restrict the numerical studies in this paper to Pythia 6, 
for completeness we also include the case of Pythia 8 | |50| , for which the corresponding relevant parameters 
are 



• TimeShower : alphaSvalue, 

• TimeShower : alphaSorder, 

• SpaceShower : alphaSvalue, 

• SpaceShower : alphaSorder, 



for final-state (timelike) and initial-state (spacelike) showers, respectively. Notice in particular that one 
here specifies the value of as{Mz) rather than that of Aqcd- Similar comments about the effective scheme 
definition as for Pythia 6 apply. 



Radiation Phase Space: The size of the allowed phase space for radiation in the shower generator may 
also affect the matched result. In the /Jx-ordered shower in Pythia 6, the switch MSTP(72) controls the 
starting scale for final-state radiation off jets that are produced by initial-state radiation and/or are colour- 
connected to the beam. Naively, the FSR off such a parton should start at the scale at which it was created, 
which is obtained with MSTP (72) =2, the recommended option. Using the other available options is strongly 
discouraged, as these lead to a quite bad agreement with NLL resummations, underscored by Banfi et al 



in 151 1. The Perugia 2011 tunes all use MSTP(72)=2. 



^It is probably closest to the so-called CMW scheme (28) . For completeness, a reasonable derived guess for a corresponding 
MS value would then be ~ 0.16, since the CMW prescription yields effective values that are approximately 1.6 times larger than 
the MS one. 
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Further, in the Perugia 2011 tunes, we start both the ISR and FSR evolutions at pTevoi = SCALUP, with the 
Pythia evolution variable pTemi defined in [4| and SCALUP the scale parameter defined in the Les Houches 
Accord for event generators |52 53|. Technically, this is achieved by setting PARP(71) = PARP(67) = 
1 . 0, where the former controls the scale factor applied to the starting scale for FSR and the latter sets the 
one for ISR. Note that these parameters could still be varied somewhat around their central values, since the 
T'revoi variable used by Pythia is not 100% identical to the Pf definition that might be used to place cuts in a 
matrix-element generator, but we have not judged this difference essential at the current level of precision. 



A.2 PYTHIA Tunes 

In this work, we have used the Perugia 0, Perugia Hard, Perugia Soft, Perugia 2010, Perugia 2011, Perugia 
2011 radHi, Perugia 2011 radLo |6|, and the DW |7| tunes of Pythia 6.4 [2]. All use the CTEQ5L PDF 



set |33|. For a complete description, see the indicated references. The salient features of the tunes are as 



follows. 

Tune DW Q is a tune of the 2^-ordered shower. It is based on the Tevatron "Tune A", which had great 
success in describing the underlying event measured at the Tevatron. Contrary to Tune A, however, DW 
also included the Drell-Yan pi spectrum, for which Tune A predicted a far too soft spectrum. Tune DW 
therefore has a significantly lower renormalization scale for ISR (and thus a larger value of a,), and 2 GeV 
of so-called "primordial k±", as compared to 1 GeV in Tune A. The energy scaling of the underlying event 
was based on comparisons between the underlying-event level at the Tevatron between 630 and 1800 GeV. 

The Perugia tunes Q are all tunes of the px-ordered shower. Unlike DW, which was developed by tuning 
to the underlying event in jet events, the Perugia tunes primarily used minimum-bias data as drivers, relying 
on the universality of PYTHIA's MPI modelling to extrapolate to the underlying event. In addition, a 
comprehensive update of the LEP fragmentation parameters was included in all tunes. 

The first set, the Perugia family, used LEP event shapes and fragmentation data, Tevatron minimum-bias 
data, and the Tevatron Drell-Yan pj spectrum. Again, the scaling from Tevatron data at 630 GeV was 
used to determine the scaling with CM energy, with some additional constraints from older UA5 data also 
included. A "Hard" and "Soft" variation attempted to vary the shower radiation up and down, respectively. 
Both Perugia and the "Hard" variation use the so-called "CMW" scheme for Aqcd for ISR, while the soft 
retained the unmodified MSbar value, in all cases taking the numerical value from the PDF set used. For the 
"Hard" variation, the renormalization scale for ISR was 0.5 pj, for Perugia pj, and for the "Soft" variation, 
V^Pj. In addition, the "Hard" variation had higher-than-nominal values for FSR, and had a slightly harder 
hadronization spectrum, while the converse was true for the "Soft" one. None of these early tunes used the 
recommended MSTP(72)=2 setting, and hence predicted rather narrow ISR jets. The Pythia tune numbers 
are 320, 321, and 322, for Perugia 0, "Soft", and "Hard", respectively. 



Tab|3]lists the parameter settings of the Perugia family that were used for the block variations in 2.2 



In Perugia 2010, jet shapes were included among the tuning constraints. The amount of FSR outside reso- 
nance decays (previously controlled by the Aqcd value read from the PDF set) was adjusted to agree with the 
level inside them (constrained by fits to LEP event shapes), combined with the recommended MSTP (72) =2. 
The Aqcd value for ISR was still read from the PDF set, and translating from MSbar to CMW, as in Perugia 
0. A few fragmentation parameters were slightly revised, since some of the previous ones had only been 
constrained using the 2^-ordered shower, and a new colour-reconnection model was introduced. No "Hard" 
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Table 3: Table of Perugia tune parameters relevant for this study. For a complete list see |[6| 

and "Soft" variations were produced for this tune. The Pythia tune number is 327 for Perugia 2010. 

In Perugia 2011, it was possible to include some early lessons from LHC at 7 TeV. Based on observed 
strangeness and baryon production rates, a few of the fragmentation parameters were again revised. The 
universal effective Aqcd choice advocated in this paper was introduced. Variations labelled "radHi" and 
"radLo" were defined as well, expressing a factor 2 variation in the Aqcd values used for ISR and FSR. The 
Pythia tune numbers for Perugia 201 1, radHi, and radLo, are 350, 35 1, and 352, respectively. 

Tabulated values of the parameters of all of the Perugia tunes can be found in the appendices of ||6|. 
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